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Abstract 

We present a detailed phcnomcnological study of forward hadron (ir°) production in deep 
inelastic scattering, with both the direct and the resolved contributions calculated to NLO 
accuracy. A comparison of the theoretical predictions for the various distributions with the 
HI data and a study of stability of the QCD predictions under changes of scales is the focus 
of this study. We obtain a very good overall description of the recent HI data with the 
choice of scale Q 2 + E\, in contrast to the (Q 2 + E^_)/2 required earlier when the resolved 
contribution was included only at LO accuracy. We find a more modest variation of the 
predictions, as the scale is changed from (Q 2 + Ej_)/2 to 2(Q 2 + E 2 ), as compared to the 
case where the resolved contribution was included only at LO accuracy. This variation is of 
the order of the rather large experimental errors. Unfortunately, this fact prevents us from 
concluding that perturbation theory gives an unambiguous prediction for forward particle 
production in deep inelastic scattering. However, the overall success of perturbative QCD in 
explaining the small XBj data means that perhaps a full resummation of the BFKL ladder 
is not called for. We notice the need for rather large resolved contributions to explain the 
data at low xbj even at somewhat larger Q 2 values. 
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1 Introduction 



Recent experimental data, from the HI collaboration pQ, on large transverse energy hadron 
production in deep-inelastic scattering have generated several theoretical papers attempting to 
explain the data within the framework of perturbative Quantum Chromodynamics in the next- 
to-leading order (NLO) approximation. These experimental results confirm and extend older 
data from HI j2j 0| and ZEUS [3J El- Since these HI and ZEUS data on forward hadrons 
as well as the data on forward jet [SJ [Tj |H1 production at large transverse momentum probe 
the small Bjorken-xej region, it was argued [Hj that they would be ideally suited to probe 
the Balitsky-Fadin-Kuraev-Lipatov ^U] (BFKL) regime, where the resummation of ln(l/xBj) 
terms is important, and that they would show the breakdown of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) regime. In this respect the single hadron data are more relevant 
than the jet data since they cover a lower Bjorken-x^j range, down to XBj = 4.10 -5 pQ. These 
data are also expected to be more accurate than the jet data because of the difficulty of jet 
identification at low transverse momentum and in the forward region. 

A comparison of older HI [2] results on single 7r° production, with a model based on lowest 
order matrix elements and parton cascades ^2], shows very strong disagreement between data 
and theory. The model falls much below the data (a factor 5 to 10 at low Q 2 and low XBj) and, 
besides, the shape of the XBj dependence is incorrect. Adding the contribution where the virtual 
photon is resolved ^3] reduces somewhat the disagreement at large Q 2 but falls short of the 
data at small Q 2 unless a very large scale is chosen in the evaluation of the anomalous photon 
component pQ. However, predictions based on improved leading order BFKL dynamics |14j 
show a better overall agreement when compared to HI data, specially at low Q 2 , but they do 
not describe the Q 2 evolution correctly pQ. 

The recent theoretical developments concern mainly the calculation of the single hadron pro- 
duction in the NLO approximation |151 1161 1171 1181 119j . Higher order diagrams neglected in the 
earlier approaches modify the picture in several ways. They generate new topologies and new 
hard scattering processes which should be considered as new Born terms. For example, at the 
lowest order (LO) the hard scattering terms are mediated by quark exchange while in the NLO 
approximation processes with gluon exchange appear and these become specially important 
in the forward region where the presence of the gluon pole enhances such terms. Also, at NLO, 
terms associated with the qq collinear component of the virtual photon, which build up the pho- 
ton structure function (the so-called resolved component), appear. As is well known, a "large" 
logarithm arises, asymptotically of type In^E^/Q 2 ) when E\ 3> Q 2 (E± is the hadron trans- 
verse momentum in the 7*-proton center of mass frame 2 ), and is associated with this structure 
function, specially when the photon virtuality is small. This term can then be considered a 
leading order term although it technically appears when calculating higher order diagrams |2L)| . 
This is the reason why it was introduced in ^H] where it indeed helped improve agreement with 
the data. However, using a large scale in the photon structure function to enhance the resolved 
contribution appears artificial. Indeed in an NLO calculation, the increase of the Born resolved 
contribution is compensated by a decrease of the higher order direct contribution, not included 
in ^3], in such a way that the sum is more stable under changes of scales. 

It should be stressed that single hadron production in deep-inelastic scattering experiments 
present a very stringent consistency test of perturbative QCD and its various input distributions: 
it involves the proton structure function as well as the hadronic fragmentation functions, all 
quantities rather precisely measured in other experiments. As just discussed it is also very 

2 The variable which we denote E± here is called in Ref. pQ. 



sensitive to the virtual photon structure function |211 1221 I23j . which is less well known but 
which has been recently discussed in detail in |19j . 

The common features of recent NLO calculations |151 1161 1171 118j are the following. The cross 
section contains two (complicated) pieces: the "direct" cross section, where the virtual photon 
couples directly to the hard process, and the "resolved" cross section, where the photon acts as 
a composite object which is a source of collinear partons taking part in the hard subprocess. 
The direct contribution is calculated in the NLO approximation, i.e up to 0(a 2 ), while the 
cross section involving the resolved component is calculated to lowest order (LO) accuracy with 
the photon scale compensating term included in the higher order part of the direct piece. No 
DGLAP type resummation is performed on the virtual photon structure function. Using modern 
proton structure functions |241 125j and fragmentation functions a very good agreement is 
achieved with the data when using a common (renormalization, factorization and fragmentation) 
scale set equal to (Q 2 + E 2 _)/2. However all the above calculations exhibit the same large 
scale dependence of the predictions mainly associated with the renormalization scale as will 
be seen below. A rather large sensitivity of the predictions to the fragmentation functions is 
also observed, with the data clearly favoring |TJj|, like other hadronic data, the parametrization 
of Kniehl, Kramer and Potter |2H] (KKP) over that of Kretzer [271 . Furthermore, in ^Fj\ a 
discussion is given to isolate the origin of the large corrections terms and it is found that they 
are associated with processes with a gluon exchange which are interpreted as the Born terms of 
the BFKL ladder. The theoretical papers differ in the procedure to obtain the cross section: in 
[151 116j a calculation of the single particle spectrum is performed with the infrared divergences 
compensated analytically while in |17| I18j a Monte-Carlo generator at the partonic level is 
constructed with a numerical compensation of divergences. 

In ^S]) the first evaluation, at the NLO accuracy, of the resolved contribution is presented: it 
includes both the construction and the use of the NLO virtual photon structure function as well 
as the NLO calculation of the hard matrix elements for the resolved processes. In the limited 
phenomenological analysis performed, good agreement with the data is obtained with the scale 
(Q 2 + Ej), larger than that of the previous NLO calculations. The importance of the resolved 
contribution to the cross section is again emphasized and it is shown that its factorization scale 
dependence is reduced at the NLO accuracy compared to the LO calculation. It is then expected 
that the full cross section will be less scale sensitive than in the work of [151 1161 IT71 118j . 

In the following we present a detailed phenomenological study of hadron production in deep- 
inelastic scattering with both the direct and resolved contributions calculated at NLO accuracy. 
A special emphasis will be put on the study of the scale variation of the cross section to determine 
the domain where the perturbative QCD approach is reliable, i.e. stable under changes of scales. 

In the next section we set up the theoretical framework and discuss the instabilities related to the 
various scales introduced in the calculation (factorization scale M on the proton side and M 7 on 
the photon side, fragmentation scale Mp and factorization scale [i). A detailed comparison with 
the various experimental HI distributions 1 is performed next: at small Bjorken-x^j the large 
corrections are found to be related to BFKL-like terms which appear in the NLO calculation in 
some approximation (first corrections in a s ln(l/x£j))- Studying the Q 2 dependence of the cross 
section will probe the photon structure function as it is expected to play a dominant role at low 
Q 2 while at large Q 2 the direct term is expected to dominate. Finally, studying the rapidity r]^ 
or x n = E l ® b /E l p b distributions, as well as the transverse momentum distribution of the pion 
will help constrain the quark and gluon fragmentation into pions. 

Hadron production in DIS experiments offers a very rich structure: it is a two scale problem, 
Q 2 and E? , with a large variation in the ratio of these two scales allowing the testing of the 



theoretical results in different regimes. Furthermore, the small Q 2 limit makes it possible to make 
contact with photoproduction experiments. Combining all the data will help in understanding 
the non-perturbative input to the photon structure function and its decreasing importance when 
the virtuality of the photon increases. It will give some insight on the transition from the non- 
perturbative to the perturbative regime, for the photon structure function. 

2 Theoretical framework 

We first discuss the features of the resolved cross section which were not taken into account in the 
previous papers ^HIEIIEIE]- When calculating the higher order (HO) corrections to the direct 
contribution there appear configurations where the virtual photon turns into an almost collinear 
q-q pair with the quark or the antiquark subsequently interacting with a parton from the proton. 
This HO contribution, in principle negligible for E\ close to Q 2 , is important when E 2 ^ becomes 
large. In this case we cannot content ourselves with the lowest order expression of the quark 
distribution in the virtual photon, proportional to ln(Ej^/Q 2 ). The latter must be replaced 
by a resummed LO or NLO expression. The standard procedure consists in subtracting from 
the HO corrections a term proportional to ln(M 2 / Q 2 ) and to calculate a resolved contribution 
with fully evolved parton distributions at the factorization scale M 2 . The M 2 -dependence of 
the resolved part is partly compensated by the l^E^/M 2 ) counterterm which remains in the 
direct HO contribution. In Ref. it was argued that a physical choice for the factorization 
scale is M 2 = (Q 2 + C 2 E 2 ^) with C 1 of order 0(1). With this scale, the resolved component is 
negligible when Ej_ <C Q 2 , whereas it is large when E\ 3> Q 2 . In the latter case, we recover the 
standard factorization scale C 2 E 2 ^ of large- E± reactions. 

Direct and resolved cross sections, calculated in the NLO approximation, have been discussed 
respectively in Ref. and ^H]- Here, we do not give the technical details of these calculations 
which can be found in the relevant references, but we summarize the main results that were 
obtained: 

1) The HO corrections to the direct cross section are very large (in the HI kinematical domain) 
and essentially come from graphs containing the exchange of a gluon in the i-channel. These 
graphs represent a zeroth order approximation to the BFKL ladder. 

2) The NLO direct cross section strongly depends on the renormalization scale \i. 

3) With the factorization scale M 2 = Q 2 + Ej_ the NLO resolved contribution is as large as the 
NLO direct one. Therefore we have access, through this contribution, to the parton distributions 
in the virtual photon. 

4) With the "natural" scale Q 2 + E\, the total cross section is in good agreement with the HI 
data, thus suggesting that a sizeable BFKL type contribution may not be necessary to explain 
the data. 

The second point above is important because it does not allow us to make stable predictions for 
the cross section and to assess the need for other contributions of the BFKL-type. Preliminary 
studies of the renormalization scale dependence have been performed separately, for the direct 
cross section and the resolved cross section |19j . Here we would like to do a more complete 
study of the scale sensitivity of the total cross section, including also the effects of the various 
factorization scales. As is well known, only the total cross section has a physical meaning. The 
separate contributions, Born terms, HO terms, direct or resolved terms are all factorization and 
renormalization scale dependent. Since the direct and resolved NLO cross section are available, 
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Figure 1: The cross section da/dxBj corresponding to the range 4.5 GeV 2 < Q 2 < 15 GeV 2 and 
E± > 2.5 GeV compared to HI data PP. Cuts on all the other kinematical variables are given 
in the text. The symbol HO s denotes the direct HO correction from which the lowest order 
resolved contribution has been subtracted. 

such a study of the scale sensitivity is now feasible. 

Before starting this study let us specify the various building blocks of the total cross section. For 
the parton distributions in the proton, we use the CTEQ6M tables |25| . and for the distributions 
in the virtual photon the parametrization given in Ref. |19j . In the latter this parametrization 
was used with fixed value of Q 2 corresponding to the average value <Q 2 > observed in a cross 
section. For instance for the cross section da/dxBj in the range 4.5 GeV 2 < Q 2 < 15 GeV 2 (see 
Fig- 0), we used <Q 2 >= 8 GeV 2 . This value corresponds to the overall bin 1.1 10~ 4 < XBj < 
11.0 10 -4 . However this value changes with XBj and the description of the whole XBj domain by 
a single value <Q 2 > is not accurate. Therefore in this paper we use a parametrization depending 
continuously on XBj, Q 2 and M 2 . We work in the MS renormalization and factorization schemes 
and all the scales are equal to (Q 2 + Ej_). We take nj = 4 flavors and for a s (n) we use an 
exact solution of the two-loop renormalization group equation with Ajjg = 326 MeV. The 
fragmentation functions of the partons in 7r° are those of Ref. With these inputs we 

obtain the cross section displayed in Fig. ^ and compared with HI data measured in the range 
4.5 GeV 2 < Q 2 < 15 GeV 2 Q]. Our calculations are performed at VS = 300.3 GeV and the 
forward-7r° cross section is defined with the following cuts. In the laboratory system a 7r° is 
observed in the forward direction with 5° < 9 n < 25° ; the laboratory momentum of the pion is 
constrained by x w = E l ® b /E l p b > 0.1, and an extra cut is put on the 7r° transverse momentum 
in the 7* — p center of mass system: E± > 2.5 GeV. The inelasticity y = Q 2 /xBjS is restricted 
to the range 0.1 < y < 0.6. 

We clearly observe in this figure the points 1), 3) and 4) mentioned above, and in particular, 
the very large HO s (the index s means that the lowest order resolved component has been 
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Figure 2: Cross section variations with Cj which has been defined in the text. The cross sections 
are normalized to 1.0 at Ci = 1. The choice of kinematic conditions is the same as in Fig.l. 

subtracted from the HO corrections to the direct term, as discussed at the beginning of this 
section). We also notice the importance of the HO corrections to the resolved cross section. 
With respect to the corresponding figure of Ref. [TJJ|, we note that the resolved component 
is larger at small xsj and smaller at large XBj, which improves agreement with data at large 
XBj- This is due to the fact that the average <Q 2 > is smaller at small XBj than at large xsj- 
Figure H is the starting point of our scale studies. We choose a kinematic region for which the 
HO corrections are large (this is due to the small value of E± > 2.5 GeV) in order to better 
exhibit the scale dependence, but with the consequence (as we shall see) that the cross section 
does not stabilize for an optimum choice of the scales. Let us define the factorization scales 
M% = C 2 (Q 2 -\-E"j_) where k stands for I (the proton distribution scale) or F (the fragmentation 
function (FF) scale). We also introduce the virtual photon factorization scale M 2 = Q 2 + C 2 E\ 
and the renormalization scale /j, 2 = C 2 (Q 2 + E'j). 

We study the sensitivity of the various components of Fig. 1 in the single bin 1.1 x 10 -4 < 
XBj < 11-0 x 1CP 4 and start with the factorization scales Mj and Mp. In Figs. 2 and 3, 
we observe very different behavior. The variation with Mj is almost flat whereas that with 
Mp is strongly decreasing. These differences are due to the different average values of x p , the 
proton distribution variable, and z, the fragmentation function variable, corresponding to the 
kinematics of Fig. 1. In the direct process for instance, we have <x p >~ 0.1, a domain in which 
the proton distribution functions do not vary much with Mj. For the fragmentation variable 
we have <z>~ 0.3 in the direct case, and <z>~ 0.7 in the resolved case. In these ranges, the 
variation of the fragmentation functions D(z,Mp) are not negligible; the higher the value of z, 
the stronger the variation. Hence the different behavior of the direct and resolved contributions. 
Moreover in the direct case, we observe that the HO corrections do not compensate the Born 
term variation. This is due to the fact that the HO corrections contain new channels which 
appear as new Born contributions for which there are no compensating ln(Mp/E±) terms. 
For instance, we have the opening of the new channel corresponding to Fig. 4a with the final 
gluon fragmenting into a 7r°. At NLO there is no counter term which corrects this new Born 
contribution. Such terms would only appear at NNLO. The contribution of this new channel, 
involving the exchange of a gluon in the £-channel, is very large and the overall behavior of the 




Figure 4: Examples of HO graphs leading to the opening of new channels when the final hadron 
is a fragment of the gluon or of the quark q. 




Figure 5: Cross section variations with C 7 which has been defined in the text. The cross sections 
are normalized to 1.0 at C 7 = 1. The choice of kinematic conditions is the same as in Fig.l. 



NLO cross section is very similar to that of the Born contribution. In the resolved contribution 
a graph with a gluon exchanged in the t-channel already exists at the Born level and the HO 
corrections contain the appropriate counter term. However the compensation between the Born 
contribution and the HO is not complete due to the large values of <z> involved. 

Let us now turn to the study of the variations with M 7 which are displayed in Fig. 5. In 
the variations studied till now, there was no compensation between the resolved and the direct 
terms. For instance, the Mj dependence of the parton distributions in the photon were separately 
compensated by \n(Mj /E±) terms which appear in the direct or in the resolved HO corrections. 
However, for M 7 , we have compensation between the resolved and the direct terms that we can 
observe in Fig. 5. When M 7 increases the Born+ HO s direct contribution decreases. This is 
due to the fact that a term proportional to log(M 7 /Q) is subtracted from the HO corrections 
leaving a piece log(Sj_/M 7 ) in the remaining HO s part, as explained at the beginning of this 
section. On the other hand, the resolved contribution increases with the increase of the parton 
distributions in the virtual photon. A counter term present in the HO resolved correction 
dampens the variation of the NLO cross section compared to the Born case. More precisely, 
the scale variation of the photon structure function contains two pieces (see i.e. eq. (16) in 
Ref. ;19j: the inhomogeneous part, proportional to a, and the homogeneous or hadron-like 
part proportional to aa s . The scale variation of the inhomogeneous part is compensated by the 
HO,, direct term, while that of the homogeneous part is in the HO resolved contribution. For 
a consistent calculation it is therefore necessary to work at the NLO level for both direct and 
resolved pieces. Due to the compensation between the direct and the resolved contributions, the 
total cross section exhibits a smoother behavior when varies by a factor 20. 

Finally let us consider the variations as a function of the renormalization scale /i. They are the 
largest. We note the same phenomenon as observed for the M^-scale variation: no compensation 
for the direct NLO term and a small compensation for the NLO resolved contribution. Con- 
cerning the direct term, this behaviour again arises due to the opening of new channels, without 
virtual corrections (they appear only at NNLO), containing terms in log(fi/E±) to compensate 
the fi dependence of a s (fj,). As these new Born terms are proportional to a?(/i) and constitute 
a large part of the cross section ^7j, the variation of the latter is strong. On the other hand 




Figure 6: Cross section variations with which has been defined in the text. The cross 
sections are normalized to one at = 1. The choice of kinematic conditions is the same as 
in Fig.l except for the rightmost panel, where the dash-dotted curve corresponds to the NLO 
cross section calculated in the bin 20 GeV 2 < Q 2 < 70 GeV 2 , 3.9 x 10 -4 < x Bj < 6.3 x 10~ 3 , 
E± > 3.5 GeV. 

virtual corrections to terms containing a gluon exchanged in the t-channel are present in the 
resolved contribution. This produces the small effect observed in Fig. 5 and we do not find 
any reasonable value of \x for which the cross section would reach an optimum. This is due to 
the large HO corrections corresponding to the small values of the transverse energy in the HI 
kinematical domain studied here. Indeed the HI experiment puts a minimum cut-off on E± (the 
transverse energy in the 7*-proton center of mass frame) of 2.5 GeV. This is a small value for a 
"large-pr" reaction and the resulting HO are large. However for higher values of the cut-off, the 
HO corrections are smaller and we find a //-variation of the resolved cross section which exhibits 
an optimum (maximum) point. For instance in Ref. JH] a cut-off E± > 5 GeV was used and an 
optimum of the cross section was found for ~ 0.2. 

Therefore we reach the conclusion that the addition of the NLO resolved component improves 
the behavior of the cross section with respect to the scale variation. However the sensitivity 
of the cross section to the renormalization scale variation prevents us from predicting absolute 
values for the latter. For instance, in the range 1/4 < C 2 < 4, the predictions vary by a factor 
2. This fact clearly points towards the necessity of calculating NNLO corrections. For the time 
being with the aim of phenomenological applications in mind, we choose scales which lead to a 
good description of the data in the range 4.5 GeV 2 < Q 2 < 15 GeV 2 . As we can see from Fig. 1, 
such an agreement is found with all scales set equal to Q 2 + Ej_. Then with the same scales 
we make predictions for da/dxsj in the other Q 2 -ranges, as well as for da/dE± and da/dx^. 
Because of the marked scale sensitivity of the cross sections, scales giving a good description 
of data in a given Q 2 -range do not necessarily lead to a good agreement in another range. It 
turns out, as we shall see, that a satisfactory description of all the data can be obtained with 
this single choice of scales. Of course the scale choice could be refined in order to improve the 
agreement between data and theory in Fig. 1, especially at small xsj • But this is a formal 
exercise that does not present any physical interest. 

Finally, to ameliorate some of our negative conclusions on the scale dependence of the cross 
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Figure 7: Inclusive 7r° cross section as a function of XBj in the range E± > 2.5 GeV for three 
different intervals in Q 2 . The cuts on other variables are given in the text. The data points are 
from the HI collaboration pQ. The histograms are the theoretical predictions obtained with all 
scales set equal to (Q 2 +E 2 J; solid line: full NLO predictions; dashed line: "direct" contribution. 

section, we note that at larger Q 2 and Ej_, the sensitivity to the renormalization scale is reduced. 
In Fig. 6 (rightmost panel), we display the behavior of the total NLO cross section in the range 
20 GeV 2 < Q 2 < 70 GeV 2 with 3.9 x 10~ 4 < x Bj < 6.3 x lO" 3 (E ± > 3.5 GeV). The cross 
section varies by less than ± 25% when C M is in the range 1/4 < C 2 < 4. 



3 Comparison to HI data 



We are now ready to compare the theoretical predictions to the HI recent results [T. on single 
7T° inclusive cross section. The same kinematical cuts as in the experimental data are imposed 
on the theory, and they are given in the previous section while discussing Fig. ^ 



Concerning the theoretical predictions "NLO" will refer to the full next-to-leading logarithmic 




Figure 8: Inclusive tt cross section as a function of XBj in the range E±_ > 3.5 GeV for three 
different intervals in Q 2 . The cuts on other variables are given in the text. The data points 
are from the HI collaboration The histograms are the theoretical results: the solid line 
corresponds to full NLO predictions and the dashed line to the "direct" contribution. Choice of 
scales is as in Fig. |7| 



predictions for the direct term as well as for the resolved term, where "direct" refers to the 
lowest order ("Born") term with the attached higher order corrections labeled "HO s " above 3 . 
The MS scheme is used throughout with A^^r = 326 MeV. For convenience we recall here the 
basic ingredients entering the calculation. All predictions are made using CTEQ6M for the 
proton parton distributions and KKP |26j for the fragmentation functions of the pion. For the 
virtual photon structure function, in the resolved term, the recent parametrization of Ref. 
is taken. Using the lowest order approximation (the so-called "box" approximation) changes the 
results by less than 10%. The common scale is chosen to be (Q 2 + E\). 

The comparison between theory and experiment for the single pion distribution as a function 

3 In |18| an extensive discussion is given of the interference terms where the photon couples to two different 
quark lines, which leads to triangle graphs when calculating the cross section (the so-called Furry terms) and these 
terms are found to give an appreciable contribution. In the present calculation, valid for neutral pion production, 
these terms are not present because the quark production cross section is cancelled by the antiquark cross section. 
The same reason makes them vanish in jet production. 



of xsj is shown in Figs. and |S] for the cuts E± > 2.5 GeV and E± > 3.5 GeV respectively. 
We notice the very good overall agreement between data and theory (note the linear scale) 
for the whole xsj range. At a finer level one may observe some systematics in Fig. [7J where 
at low XBj, for the medium Q 2 range, the theoretical predictions fall slightly above the data 
while in the large Q 2 bin it is the opposite. Furthermore, one notes the importance of the 
resolved contribution (the difference between the solid and the dashed line). At low Q 2 (upper 
panels) it is 1.5 to 1.9 times the direct contribution, decreasing as XBj increases, while at large 
Q 2 (lower panels) it never exceeds the direct term and becomes almost negligible at large xsy 
this is as expected from the discussion in the previous section. The importance of the resolved 
contribution to obtain agreement with the data was also pointed out by Kramer and Potter who 
calculated the NLO cross section (resolved term at leading order) to forward dijet production 28 
and compared it with HI 3 a and ZEUS [7J data, as well by Jung and collaborators 29 a in their 
analysis using a LO calculation. One may comment again on the rather unusual situation at 
low XBj where the HO s correction to the direct term can be up to an order of magnitude larger 
than the Born term (see Fig. ^) due to the appearance of the BFKL-like terms of Fig. 0] with 
gluon poles. At large XBj however one recovers the "usual" situation where the HO s piece is of 
the same order of magnitude as the lowest order term. 

A very impressive agreement is also achieved, in Fig. |UJ for the E± spectrum for all Q 2 values. 
The resolved contribution decreases with Q 2 but it remains important for all values of Q 2 and all 
transverse momenta. The E±_ distribution should be sensitive to the choice of the fragmentation 
functions and it is interesting to try different sets, in particular that of Kretzer [27]. We do not 
do it here as it has already been shown by Daleo et al. 16_ that the parametrization of (27J 
leads to predictions which fall below the data. This confirms previous studies [50 showing that 
the fragmentation functions of (27J systematically underestimate particle production in hadronic 
reactions. 

Similarly, the longitudinal momentum distribution of the pion is in remarkable agreement with 
the data both for specific Q 2 bins (Fig. ITU]) or specific XBj bins (Fig. ITTj) . Again the resolved 
component is important over the whole x n range but, clearly, it gives a decreasing contribution 
as Q 2 increases. One has to note that there is little correlation between x w and the fragmentation 
variable z: for instance, in the resolved case <z>~ .5 (calculated with the Born term only in 
the region 2 GeV 2 < Q 2 < 4.5 GeV 2 ) varies by less than 10% when x w varies between the first 
and the last bin of Fig. 1101 Therefore we cannot rely on the x n spectrum to constrain the z 
shape of the fragmentation functions. 

From the comparison with data we can conclude that perturbative QCD, in the NLO approxi- 
mation, gives unexpectedly good results, especially in view of the initial discrepancy observed at 
leading order between theory and data. Two ingredients explain this fact: the unusually large 
correction to the direct term, specially at low XBj, and the importance of the resolved photon 
contribution including the associated higher order corrections. Unfortunately, none of the in- 
clusive observables discussed here allows for an unambiguous separation of the two terms. In 
principle, looking at more exclusive quantities, such as hadron-jet correlations, would allow the 
determination of the longitudinal momentum fraction in the photon x^ |5J , and consequently 
the separation of the two types of terms. However, since the HO s term is very large, it may lead 
to a large contribution at x 7 ^ 1 making the separation from the resolved term difficult. 

The success of perturbative QCD to explain the data at small XBj is interesting. It seems to 
imply that there is no clear signal in the HI data of the BFKL type resummation effects and 
that keeping only the lowest order term in the usual perturbative sense is justified. One reason 
may be the following. The BFKL result is derived for asymptotic energies. However, at HERA 
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Figure 9: Inclusive 7r° cross section as a function of the 7r° transverse momentum in the 7*- 
proton center of mass frame for three different intervals in Q 2 . The data points are from the 
HI collaboration pQ. The histograms are the theoretical results: the solid line corresponds to 
full NLO predictions and the dashed line to the "direct" contribution. Choice of scales is as in 
Fig. The variable in the figure is the notation of the HI collaboration and is called E± in 
the text. 



the rapidity range, hi{S/Q 2 ), available is not extremely large and threshold effects do not allow 
for the full formation of the BFKL ladder [31133]. 



4 Comparison with other perturbative calculations 

The HI data are also in very good agreement with the NLO calculations of Daleo et al. ^0] and 
Kniehl et al. ^H] • We recall that the difference between these approaches and the present one lies 
in the fact that, in the former, no special consideration is given to the photon structure function: 
the NLO correction to the direct term contains a "large" factor of type ln((Q 2 + E 2 J/Q 2 ) which 
amounts, in fact, to parametrizing the photon function by its lowest perturbative approximation. 
Furthermore no NLO corrections are included in the resolved cross section. In contrast, in this 
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Figure 10: Inclusive it cross section as a function of x w = E l ® /E^ in the range E± > 2.5 
GeV for three different intervals in Q 2 . The cuts on other variables are given in the text. The 
data points are from the HI collaboration pp. The histograms are the theoretical results:the 
solid line corresponds to full NLO predictions and the dashed line to the "direct" contribution. 
Choice of scales is as in Fig. 

work, we use both a NLO expression for the resolved photon structure function and we include 
the HO corrections to the resolved cross section. Agreement with the data is obtained in all 
cases at the cost of using a different choice for the common scale. In ^S] and ^S], as well as 
in our previous work the scale (Q 2 + Ej_)/2 was the appropriate choice. In this work, it is 
seen that the scale (Q 2 + Ej_ ) is preferred. The data do not obviously prefer one or the other of 
the two sets of calculations as the shape of the observables is not affected. In Ref. ^H] a very 
large scale sensitivity was however observed: under the rather modest change from (Q 2 + £' 2 _)/4 
to (Q 2 + E"j_) the theoretical predictions vary by as much as a factor 2 in some cases, and, in 
any case, the theoretical uncertainties are much (sometimes twice) larger that the experimental 
ones (statistic and systematic errors combined). In the current approach we expect a smaller 
sensitivity to the scales since more HO corrections are taken into account. Besides, it is seen 
from Figs. |2J EJ El that the variation with the scales seems to decrease at higher scales. This is 
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Figure 11: Inclusive it cross section as a function of = E l ® /E ^ in the range E± > 2.5 
GeV for three different intervals in xsj- The cuts on other variables are given in the text. The 
data points are from the HI collaboration 1 . The histograms are the theoretical results: the 
solid line corresponds to full NLO predictions and the dashed line to the "direct" contribution. 
Choice of scales is as in Fig. |7| 

illustrated in Figs. WI\ and IT31 where we show the results for (Q 2 + E\)/2 and 2(Q 2 + E\) 4 . 
Compared to the results of JH]) the scale variations are somewhat tempered and are of the same 
order as that of the rather large experimental errors. 



5 Conclusions 

Using the latest structure and fragmentation functions, the complete NLO calculation of the 
direct and resolved contributions to forward particle production in deep-inelastic scattering at 
HERA, describes the data rather well in the wide kinematical range available: 2 GeV 2 < Q 2 < 

4 Note that the resolved photon scale differs slightly from the others since we use, as explained in Sec. [5] 
Af2 = (Q 2 + C 2 E\) with C 2 = 1/2 and C 2 = 2. 




Figure 12: Inclusive 7r cross section as a function of XBj in the range E±_ > 2.5 GeV for three 
different intervals in Q 2 . The data points are from the HI collaboration pQ. The histograms are 
the NLO theoretical results for different scales of the form C 2 (Q 2 + E 2 ^): C 2 = .5, upper dotted 
histogram; C 2 = 1, solid histogram; C 2 = 2 lower dashed histogram. 



70 GeV 2 , 2.5 GeV < E± < 15 GeV. The importance of the NLO corrections to both the direct 
and resolved terms is pointed out. These large corrections are associated with new topologies 
involving gluon exchange in the hard sub-processes. These terms, which have no equivalent 
at the lowest order are interpreted as the first terms of the BFKL ladder. The data seem to 
indicate that resummation of such ladder diagrams is not necessary, probably because of the not 
so large rapidity phase space available. Agreement between theory and data is achieved choosing 
a standard scale of the form (Q 2 + E 2 ). The variations under the proton factorization scale 
and the photon factorization scale are under control. However a rather large instability of the 
predictions is observed when varying independently the renormalization and the fragmentation 
scales. This prevents a really quantitative prediction for the single pion inclusive distribution in 
the forward region. In this respect, taking account of the HO resolved contribution improves the 
situation compared to calculations which ignored it but the situation is still far from satisfactory. 
We have checked that imposing a larger E±_ cut on the data reduces the scale sensitivity: for 
example, with the conditions E± > 7 GeV and 4.5 GeV 2 < Q 2 < 15 GeV 2 the variation is ±^ 
for a scale variation as in Fig. 1121 Probably, the evaluation of the next-to-next-to-leading order 
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Figure 13: Inclusive ir° cross section as a function of xsj in the range E± > 3.5 GeV for three 
different intervals in Q 2 . The data points are from the HI collaboration ^p. The histograms are 
the NLO theoretical results for different scales of the form C 2 (Q 2 + E 2 ): C 2 = .5, upper dotted 
histogram; C 2 = 1, solid histogram; C 2 = 2 lower dashed histogram. 



terms is required to obtain reliable and stable perturbative predictions. 
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